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ABSTRACT 

We investigate the structure and stability of hypercritical accretion flows 
around stellar-mass black holes, taking into account neutrino cooling, lepton con- 
servation, and firstly a realistic equation of state in order to properly treat the 
dissociation of nuclei. We obtain the radial distributions of physical properties, 
such as density, temperature and electron fraction, for various mass accretion 
rates 0.1 ~ lOM s _1 . We find that, depending on mass accretion rates, different 
physics affect considerably the structure of the disk; most important physics is (1) 
the photodissociation of nuclei around r ~ 100r 5 for relatively low mass accretion 
rates (M ~ 0.01 — O.1M s _1 ), (2) efficient neutrino cooling around r ~ 10— 100r g 
for moderately high mass accretion rate (M ~ 0.2 — 1.0M Q s _1 ), and (3) neutrino 
trapping (r ~ 3 — 10r 9 ) for very high mass accretion rate (M > 2.0M o s~ 1 ). 
We also investigate the stability of hypercritical accretion flows by drawing the 
thermal equilibrium curves, and find that efficient neutrino cooling makes the 
accretion flows rather stable against both thermal and viscous modes. 

Subject headings: accretion, accretion disks — black hole physics — gamma-ray 
bursts — neutrinos — instability 



1. INTRODUCTION 



Gamma-Ray Bursts (GRBs) are the most energetic explosions in Universe, which are be- 
lieved to release energies up to 10 51_53 ergs in only a few tens of seconds. Their observational 
behavior (lightcurves of prompt emission and afterglow, spectrum, etc.) are well explained 
by the relativistic fireball shock model (for reviews, see e.g. Piran 1999, and Meszaros 2002). 
The central engine which produces hot and baryon-poor plasma (i.e. a fireball) is, how- 
ever, totally hidden from direct observation, and so the engine and fireball-making processes 
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have not well been understood yet. Nevertheless, there is one popular model for the en- 
ergy sources of GRBs which has been intensely studied by now — a hyperaccreting black hole 
model (Narayan et al. 1992, 2001). In this model, the energy of relativistic jets that make 
intense gamma ray emissions is produced via the accretion of a massive disk (~ O.1M ) 
onto a stellar mass black hole. Such systems are expected after several energetic phenom- 
ena; mergers of double neutron star binaries (Eichler et al. 1989), neutron star/black hole 
binaries (Paczyhski 1991), white dwarf/black hole binaries (Fryer et al. 1999), black hole/He 
star binaries (Fryer and Woosley 1998) and failed supernovae (or collapsars, Woosley 1993, 
Paczyhski 1998, MacFadyen and Woosley 1999). 

There are several reasons why the hyperaccreting black hole model is widely accepted. 
First, this model can explain both the variability timescale (5t ~ 1msec) and the dura- 
tion timescale (T dur ~ 1 — lOOsec); the former timescale represents the compactness of the 
central engine-the inner radius of the disk is smaller than c5t ~ 3 x 10 7 cm — and the lat- 
ter represents the accretion timescale — the outer radius of the disk r out can be estimated by 
2~dur ~ r ont / (ac s (H I Tout)) = 1/ (cxQ(H I V out ) 2 ) , where a, H and c s are the viscosity parameter 
(see below), scale height of the disk and sound velocity, respectively. Second, the relativistic 
jets which are believed to be produced in GRBs formed in accretion disk systems, as can be 
seen in active galactic nuclei or microquasars. 

In such a massive accretion disk, photons are almost trapped because the flow become 
highly optically thick, and so radiative cooling is not efficient. However, as the flow become 
dense (~ 10 7 g cm -3 ) and hot (~ 10 10 K), thermal neutrino emissions will set in, and be- 
come the most efficient cooling process. Due to this nature, such disks (or flows) are called 
"neutrino-dominated accretion flows " (NDAFs). The NDAF model can give explanations for 
not only the energetics of GRBs but also the processes of making the relativistic and baryon- 
poor fireballs through the energy release by neutrino annihilation above the disk (Popham et 
al. 1999). In order to predict the amount of energy deposited through this process precisely, 
it is necessary to solve the disk structure taking into account various kinds of microphysics 
because when calculating total luminosity of neutrino annihilation we should know where 
and how much neutrinos are emitted from the disk. First, in the region with a moderate 
density and temperature, the fraction of heavy nuclei (e.g. 56 Fe, 62 Ni, etc.) can be high. 
However, around r ~ 100r 9 these nuclei are photodissociated because the temperature of the 
disk rises. This process consumes the internal energy of the accretion flow and so effectively 
works as cooling. Secondary, in such a high density disk electrons are highly degenerated 
and this will affect the equation of state (EoS), because electron degeneracy pressure will 
contribute to total pressure, and neutrino emission, because degenerate electrons tend to 
be captured by protons and electron neutrino emissivity becomes larger. Especially, by the 
latter effect the accretion flow are supposed to be highly neutron-rich. Finally, the accretion 
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flow can be optically thick even for neutrinos in the innermost region, so the energy and 
number density of neutrinos will be advected with accreting matter. All of these effects 
should be taken into account when we analyze the structure of the disk and its stability. 

There are many previous studies on NDAFs. Popham et al. (1999) performed relativis- 
ts calculations for the disk structure in Kerr spacetime with neutrino-thin approximation. 
Narayan et al. (2001) extended their discussion to the low luminosity regime in the context 
of CDAF (convection-dominated accretion flow), and Kohri & Mineshige (2002) studied the 
effects of highly degenerated electrons and nucleons on neutrino emissions and the pressure 
in the disk. Di Matteo et al. (2002) and Kohri et al. (2005) considered the neutrino trap- 
ping effect, which is not negligible when mass accretion rate is very high. Gu et al. (2006) 
discussed the contribution from the neutrino-thick region to neutrino emissions and showed 
that sufficiently large amount of energy for GRBs would be released by neutrino pair an- 
nihilation. Recently, Chen & Beloborodov (2006) calculated the structure of hyper-critical 
accretion flows around Kerr black holes, fully taking into account the effect of general relativ- 
ity. In their calculations, however, the detailed physics shown above are not fully considered 
but only part of them approximately. In this paper, we take into account the effects of 
neutrino emission and trapping in NDAFs, and employ the realistic EOS in order to see the 
nuclear composition and dissociation in the outer region of the flow consistently. 

The plan of this paper is as follows. In Section 2, we show the fundamental equations 
and assumptions describing the stational behavior of NDAFs. In Section 3, we show the 
calculational results and implications of them. Finally in Section 4. we discuss our results 
and compare them with previous studies. 

2. MODELS 

In this section, we present the basic equations and assumptions that describe the situ- 
ations of our interest. We assume that the accretion flow we are considering is steady and 
axisymmetric. In addition, we assume that the gravity is described with pseudo-Newtonian 
potential (Paczyhski and Wiita, 1980), i.e. 



GM 
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and the angular velocity and specific angular momentum of gas particles in circular motion 
are assumed to be Keplerian, i.e. 



GM 

l(r) = r 2 Q{r), (3) 

The expressions for mass conservation, angular momentum conservation, energy conserva- 
tion, lepton fraction conservation, hydrostatic balance and a viscosity are, respectively, 
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where r, E, /(r), p, H, v r , u, Q, p, e and Yi denote the radius, surface mass density, specific 
angular momentum at the radius r, baryon mass density, scale height, radial velocity, kine- 
matic viscosity coefficient, angular velocity, pressure, internal energy per baryon, and lepton 
fraction, respectively; M, a, m p are the mass accretion rate, viscosity parameter and mass 
of a proton, respectively; Q + , Q~ , F n ^ e and F n ^ e are the viscous heating rate per unit area, 
cooling rate per unit area, neutrino number flux, anti-neutrino number flux, respectively. We 
set the inner radius of the accretion disk, r in , as the radius of the innermost stable circular 
orbit of the central black hole; r in = 3r g . The surface mass density is defined as 

E = 2pH, (10) 

and the dissipated energy per unit area is 

r-*(r«P) 2 . (ID 



dr 

using the above expression of Q(r). On the other hand, the cooling rate per unit area is 
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where Q v . is the cooling rate per unit area via i-neutrino emission, each of which is calculated 
below. 

In calculating the neutrino energy flux and the neutrino number flux, we take into ac- 
count the various kinds of reaction; electron-positron capture onto nucleons/nuclei, electron- 
positron pair annihilation, and plasmon decay. The emissivities and reaction rates of these 
processes are calculated based on the formulae given by Bruenn (1985), Ruffert et al. (1996), 
and Burrows & Thompson (2002). In order to solve the energy equation, the lepton number 
equation, and the hydrostatic equilibrium equation properly, we should know the neutrino 
energy flux and number flux on the disk surface and the neutrino energy density and number 
density in the disk. We evaluate these physical values of neutrinos by adopting the approx- 
imation by Popham & Narayan (1995) and Di Matteo et al. (2002). According to their 
procedure, the neutrino energy flux, energy density, number flux and number density are 
given by 

n- = 2 (c/4K,o 

(3/4)^/2 + 1/^3 + 1/(3^)]' 

u = u Vtfi (3r E>Vi + 2V3) 
3t e ^ + 2V3 + 2/t e ^ 

F = (c/4K "° 

(3/4)^/2 + 1/^3 + 1/(3^)]' 

= ^, (3r +2^) 

3r n , Vi + 2V3 + 2/r n 

respectively. Here, TE, a ,ui and T n ^ Vi are the spectral-averaged neutrino-absorption optical 
depth for energy transport and number transport, respectively, and TE, Vi and r n jV . are the to- 
tal neutrino optical depth (including scattering) for energy transport and number transport, 
respectively. In the above expressions, u Vi $ and n Uit0 are defined as 

k 4 T 4 

Uu -° = 2^¥^ F3 ^ /kBT ^ (17) 

where F m (r]) is the standard Fermi-Dirac integral given by 

dx ^rr (19) 

and ji u . represents the effective neutrino chemical potentials. The calculation procedures of 
neutrino opacities and effective neutrino chemical potentials are given in Appendix A. 
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Most of the previous studies on the structures of NDAFs, a simplified equation of state 
and nuclear compositions (i.e. all nuclei which exist in the flow are only a particles) were 
assumed and so the estimation of the photodissociation-cooling effects and the value of 
electron fraction Y e might not be accurate in the computational results. In this study, we 
employ a realistic equation of state by Lattimer & Swesty (1991), which is widely used in 
supernova simulations. By using this, we successfully take into account the effects of heavy 
nuclei-dissociation in the calculation. 

Finally, we mention the outer boundary conditions that we impose in computing the 
structures of NDAF. We assume that the inflowing gas is composed primarily of neutron-rich 
iron group nuclei (i.e. 56 Fe, 58 Fe, 62 Ni and so on), and the electron fraction Y e is ~ 0.42. 
Such compositions are often realized in the detailed calculations of the core of a massive 
star in its final stage of evolution. In addition, we assume that in the outer boundary the 
dissipated energy in the flow is totally advected: 

rripp dp\ 
p 2 dr J 



(20) 

where s is the specific entropy and £ is the parameter of the order of unity (see Chap. 8 of 
Kato et al. 1998). The exact value of £ is not needed in this calculation. 

With these boundary conditions, we solve the set of equations from the outer boundary, 
using the fourth-order Runge-Kutta method. 

3. RESULTS 

In this section we present our results of the structure of hypercritical accretion flow with 
various mass accretion rates, and of the equilibrium curves. The flow structures with typical 
parameters are shown in Fig 1-3 for relatively low (M = 0.01 — O.O8M s -1 ), moderately 
high (M = 0.1 — O.8M s _1 ), and very high (M = 1.0 — 10.0M Q s -1 ) mass accretion rates, 
respectively. We can see from them that the profiles of density and temperature are highly 
dependent on mass accretion rates. In the following, we see dominant physics in each of 
accretion regimes 
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3.1. Cases of M = 0.01 - O.O8M s" 1 

First, we show the structures of the disks with M = 0.01 — 0.08M Q s -1 in Fig 1. When 
disk temperature reaches ~ MeV around r ~ 100r 9 the photodissociation of nuclei occurs. 
We can see that the density profiles get slightly steeper, and the temperature profiles are flat- 
ter at r < 100r 9 than in the other region (Fig 1). These kinds of behavior can be understood 
as the effects of efficient cooling. In the region of photodissociation, viscously dissipated en- 
ergy is efficiently consumed rather than advected inward, and the flows are effectively cooled. 
Therefore, the disks become thinner due to lower temperature and pressure, and densities 
get higher compared with the case in which dissociation is not taken into account. The 
regions of photodissociation appear around r ~ 100r 9 with arbitrary mass accretion rates. 
In the photodissociation region, the degeneracy parameter of electrons r] e rises because of 
the photodissociation cooling. 

3.2. Cases of M = 0.1 - O.8M s" 1 

Next we show in Fig 2. the disk structures of the disk with moderately high mass 
accretion rate (M = 0.1 — O.8M s _1 ). Around r ~ 30r s , the density profiles get steeper and 
the temperature profiles get flatter, as seen in the photodissociation region. Both of these are 
due to efficient neutrino cooling. In this region, the gradient of entropy profile get positive 
(see lower-middle panel of Fig 2). Almost the same argument as in the photodissociation 
region can be applied in this region. Namely, there should be effective cooling operating. 

In Fig 4, the ratios of neutrino cooling flux (F_)to viscous heating rate per unit disk 
surface area (F + )are shown. We can see that, with moderately high mass accretion rate, 
there appears the region where F_/F + gets drastically high or even exceeds unity (i.e. larger 
amount of energy than that viscously dissipated is emitted as neutrino per unit surface area). 
Around that region, due to the sufficiently high density and temperature, thermal neutrino 
emissions become very efficient. Therefore, the accretion flow is well cooled there, and the 
profiles of density and temperature show the similar behavior as that in photodissociation 
region. Note that we have always Q ra ,d/Qvis ~ 1 in the standard disks. In NDAF, excess 
cooling is compensated by advective heating. 

Moreover, electron fraction Y e dramatically decreases inward in this region. This is 
because electron capture process (e + p — > n + v e ) becomes efficient with sufficiently high 
temperature (> lMeV) and high electron degeneracy (i] e > 1). Around there, Y e drops to 
as low as ~ 0.1, which means that the inner region of NDAF with moderate mass accretion 
rate is neutron-rich. 
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3.3. Cases of M — 1.0 - 1O.OM s" 1 

For very high mass accretion rate (> 3.OM s _1 ), we can see another interesting feature 
in the innermost region of an accretion flow, where the density profile becomes flatter, and 
the temperature profile becomes steeper than in the case with intermediate mass accretion 
rate (i.e. M ~ 0.3 — l.OM s). This means that neutrino cooling which is efficient in lower 
accretion rate disk ceases to be powerful because the flow becomes thick for neutrinos. The 
upper panel of Fig. 5 shows the profiles of the optical depth of u e in the disks with very high 
mass accretion rates. According to this figure, in the innermost region, the accretion flows 
can be thick (r„ e > l)with respect to neutrinos when mass accretion rates are larger than 
~ 3.0M Q s -1 . This makes an obvious change in the electron fraction profiles. If neutrinos are 
trapped inside the flow, neutronization (p + e~ — > n + u e ) are effectively prevented because of 
the inverse reaction. As a result, Y e , which decreases inward in the efficient neutrino-cooled 
region, no longer decreases but rather increases inward in the neutrino-trapping region, and 
takes the value around ~ 0.2 in the vicinity of inner edge. The steepness of density profiles 
or the flatness of temperature profiles are not so conspicuous as the cases with lower mass 
accretion rates. The reason of this difference can be seen in the plots of free nucleon fraction 
(the lower-right panel of Fig 3). For high mass accretion rates, the density gets high and so 
photodissociation of nuclei occurs rather gradually inward. As a result, the peculiar features 
of density and temperature profiles are smeared out. 

We can see other interesting features in the flow with even higher mass accretion rates. 
In the innermost region, entropy per baryon and electron degeneracy parameter take approx- 
imately the same values regardless of mass accretion rates. Such trends of physical quantities 
appear when the accretion flow gets thick with respect to neutrinos (see also Fig. 5). 

Before ending this subsection, we should mention the oscillations which appear con- 
spicuously on the curves of temperature, degeneracy parameter, and neutrino optical depth. 
These features can be explained in the following way. 

In the outer region of the disk, most baryons are in the form of nuclei, and the main 
opacity for the neutrinos is from the coherent scattering by nuclei. Around the radius where 
many oscillations on the curves appear, the photodissociation of nuclei occurs and the optical 
depth with respect to neutrinos drops inward suddenly. As far as the mass accretion rate is 
not so high, the optical depth is much smaller than unity even if the plenty of heavy nuclei 
exist. However, is the mass accretion rate is sufficiently high, even in the outer region the 
optical depth can be as large as the order of unity, and due to the dissociation the accretion 
flow become suddenly neutrino-thin inward. Then the neutrino cooling become efficient and 
the temperature can get lower inward, and the fraction of nuclei, which is calculated from 
EoS table, gets suddenly increased inward. This causes many oscillations on the curves of 
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some physical quantities. Note that the physical quantities in the EoS table are calculated 
assuming nuclear statistical equilibrium (NSE). In our case, this assumption is sometimes 
violated; if the timescale for nucleons to attain NSE is longer than the accretion timescale 
the composition of nuclei in the flow is expected to deviate from that evaluated from the 
EoS table. In order to calculate the nuclear composition more precisely, it is necessary to 
solve the nuclear reaction equations, which is beyond the scope of this work. 

3.4. Nuclear Composition of Hypercritical Accretion Flows 

Now that we solve the structures of hypercritical flows with arbitrary mass accretion 
rates, we can see how the composition in the disk depends on the radius and the mass 
accretion rate. In the previous studies, only a-particles and their dissociation were considered 
in the outer region of the disk, and the amount of a-particles was estimated assuming nuclear 
statistic equilibrium. This procedure, however, missed the existence of heavier nuclei (such 
as Fe, Ni, etc.), and so they could not take into account photodissociation cooling correctly, 
and could not follow the nuclear composition and the electron fraction consistently. In our 
calculations, these points are greatly improved by employing the realistic EOS by Lattimer 
& Swesty (1991)and properly considering the existence of nuclei other than He. 

Fig. 6 show the contour plots of electron fraction, a-particle fraction, and heavy nucleon 
fraction on (r, M). We can see that the existence of heavy nuclei is not negligible around 
r ~ 100r 9 . Heavy nuclei also affect the opacity of NDAFs with very high mass accretion 
rate. In the lower panel of Fig. 5, scattering optical depths with respect of electron neutrinos 
are shown with various mass accretion rates, and we can see that the scattering opacity by 
nuclei is not negligible in NDAFs with sufficiently high mass accretion rate. This fact has 
not been mentioned in previous studies, in which the existence of heavy nuclei is neglected. 

3.5. Stability Analysis of Hypercritical Accretion Flows 

We can draw the thermal equilibrium curves of NDAFs on (E, T) and (E, M) plane by 
fixing some typical radii. These curves are helpful in checking qualitively the thermal and 
viscous stability of accretion flows (e.g. Di Matteo et al. 2002). If there is an S-shaped 
sequence in the equilibrium curves on (E, T)/(E, M), the disk is thermally/viscously un- 
stable (Kato et al. 1998). We expect some features in the equilibrium curves when dominant 
physics alternates others. For example, S-shape appears when disks are cool enough for hy- 
drogens to recombine in the case of dwarf-novae. Another S-shape appears when radiation 
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pressure dominates over gas pressure. We show those equilibrium curves in Fig 7 and 8 for 
r = 10r 9 , 30r 9 , and 100r g . These figures show that the NDAFs whose structures we calculate 
are thermally and viscously stable. This is because dissociation of nuclei does not induce 
drastic changes of physics dominating the disk structure. However, we can find a remarkable 
features in the equilibrium curves. At smaller radii, the equilibrium curves on the (E, T) 
plane have negative gradient part in a certain surface mass density range, and those on the 
(E, M) plane become comparatively flat with the same E. This behavior corresponds to 
the efficient neutrino cooling in the disk with moderate mass accretion rate. As this is the 
opposite tendency compared with the case of unstable disk, this feature implies that efficient 
neutrino cooling makes NDAF rather stable both thermally and viscously. Such feature has 
not been found in the equilibrium curves of accretion flows with lower mass accretion rate 
(M~M Edd ). 

3.6. Neutrino Luminosity from Hypercritical Accretion Flows 

Fig. 9 shows the sum of neutrino luminosities of all flavors, electron neutrino luminosity 
and anti-electron neutrino luminosity (upper panel), and efficiency of neutrino emission 
(L u /Mc 2 , lower panel)as the function of mass accretion rate. According to this plot, the 
energy conversion from the accretion flow to neutrino emissions is most efficient when mass 
accretion rate is ~ O.3M sec _1 . Below this rate density and temperature in the accretion 
flow are not high enough to emit neutrinos efficiently, and above this rate the inner region 
of the disk gets thick with respect of neutrinos and so the number of neutrinos which escape 
the disk becomes smaller. 

Note that, though we have taken into account the effects of heavy nuclei in the analyses, 
the total neutrino luminosities do not differ significantly from those calculated with the same 
mass accretion rate in the previous studies, because most of neutrinos are emitted from the 
innermost region of the disk (r < 10r ff ), while the photodissociation of heavy nuclei occurs i 
the outer part of the disk (r ~ 100r g ). 

4. DISCUSSIONS 

In this paper, we investigate the structure of hypercritical accretion flows as the central 
engine of GRBs with various mass accretion rates, and the thermal and viscous stability by 
drawing the thermal equilibrium curves on (E,T) and (E, M) plane, respectively. We have 
found that the photodissociation of nuclei, efficient neutrino cooling, and neutrino trapping 
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cause some remarkable changes to the disk structure. In the calculation, we treat appro- 
priately the composition of nuclei and the neutrino optical depths by adopting a realistic 
equation of state and various kinds of neutrino reaction. By doing this, we firstly include the 
neutrino emissions and scattering by heavy nuclei, whose amount is not negligible around 
r ~ 100r 9 assuming nuclear statistical equilibrium. For relatively low mass accretion rates, 
remarkable differences from previous works do not appear in our calculations. However, for 
very high mass accretion rate, even in the outer region (r ~ 100r 9 ) neutrino optical depths 
become of the order of unity. This is because of the large opacity due to the coherent neu- 
trino scattering from heavy nuclei (See the lower panel of Fig. 5). From these facts, we 
can conclude that the existence and dissociation of heavy nuclei remarkably affect the disk 
structure and neutrino emission, and this has not been mentioned in the previous works. 

We have also investigated the neutron-to-proton ratio in hypercritical accretion flows 
by taking into account the neutrino trapped inside the flows. For moderately high mass 
accretion rates (M ~ 0.3— 1.0M Q s -1 ), because of high density and efficient photodissociation 
cooling, the electron degeneracy parameter drastically increases inward around r ~ 100r 9 . 
In addition, the temperature becomes so high that many electrons have enough energy to be 
captured by protons. Then in the inner region electron capture proceeds as the most efficient 
neutrino emission process, and then the electron fraction Y e decreases inward down to~ 0.1. 
This tendency is positive for the acceleration of an outflow that launches from around the 
innermost region of the disk; if the outflow picks up the disk material and becomes neutron- 
rich, the outflow will be accelerated more efficiently compared to the case that it picks up 
the same amount of disk material but neutron-poor, because neutrons do not interact with 
other particles electrostatically and will decouple while protons and electrons are accelerated 
(e.g. Vlahakis et al. 2003). Moreover, the disk outflow is considered to be a potential site 
for nucleosynthesis, and the conditions with large neutron-to-proton ratio is favorable to 
r-process (Surman et al. 2006). 

On the other hand, if mass accretion rates are very high (M > 3.OM s" 1 ), in the 
innermost region Y e increases inward up to ~ 0.2. This is because electron-neutrino optical 
depth becomes very high and electron neutrinos are trapped inside the accretion flow. The 
same kind of features are seen in calculations of core-collapse supernova explosion, in which 
neutrino trapping due to the coherent neutrino scattering opacity by heavy nuclei stops 
the neutronization of stellar material because neutrinos become degenerated. Moreover, as 
seen in Fig. 8, the efficiency of neutrino emission from a whole disk surface peaks around 
M ~ 0.3M Q s _1 and for larger mass accretion rates, L u /(Mc 2 ) is suppressed because of 
neutrino-thickness of the accretion flows. These results imply that the accretion disks with 
too high mass accretion rate are not always favorable for the production of outflows with 
large Lorentz factor, which are considered to be necessary to explain observed GRB emission 
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features. 

Our study also gives some implications to the possibility of energy deposition via neu- 
trino pair annihilation (v + v — > e~ + e + ) above the neutrino-cooled disk, which make a 
hot and baryon-poor fireball. In order to estimate the energy deposition rate, we should 
calculate both the energy flux of v e and v e from the disk. We investigate the effects of 
neutron-to-proton ratio, electron degeneracy, and neutrino trapping on the flavor of neutri- 
nos mainly emitted from the disk surface. As seen in Fig. 8, the luminosity of neutrinos 
from the whole disk surface is around L u ~ 10 52 erg s -1 , and there is little difference between 
the luminosities of v e and z/ e , which had often been assumed in the previous studies on neu- 
trino pair annihilations above the neutrino-dominated accretion disk (e.g. Gu et al. 2006). 
So we can say that their estimates for energy deposition rate above the disk are essentially 
unchanged. For the precise calculation, however, we should know the energy spectrum of 
emitted neutrinos because the cross section of pair annihilation depends on the energy of 
neutrinos and antineutrino. We are now investigating the possibility of neutrino spectra from 
neutrino-thick accretion disks to deform from Fermi-Dirac distribution and its influence on 
the energy deposition rates (Kawabata et al. in preparation). 

Finally, we mention the stability of hypercritical accretion flows. In this paper we have 
obtained the thermal equilibrium curves of the accretion flow at several radii where the 
dominant physics operating the flow change, and have found that hypercritical accretion 
flows are both thermal and viscously stable, even at the radius where the photodissociation 
of heavy nuclei occurs and the nuclear composition and opacity sources drastically change. 
Therefore, we cannot explain the time-variable behavior in the context of these kinds of 
instabilities. We have not considered here the effect of gravitational instability of the disk 
(Toomre 1964). As already indicated by many authors, the accretion flows with M > 
1M Q s _1 are gravitationally unstable in their outer part (e.g. Di Matteo et al. 2002). 
Such disks can fragment into gravitationally bound blobs, and this can be responsible for 
X-ray flares observed in the early phase of GRB afterglows (Perna et al. 2006). In order to 
investigate the behavior of fragmented hypercritical accretion flows, one should reconstruct 
the time-dependent disk model taking into account the self-gravity, which is left for future 
work. 
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We also thank S. Nagataki, S. Uematsu, and R. Takahashi for useful discussions and help- 
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A. Computations of Neutrino Emission Rates and Opacities 

In this appendix we give the way to calculate the neutrino emission rates and neutrino 
optical depth in an NDAF along the vertical direction. Most of the expressions we show 
below are covered by Bruenn (1985), Ruffert et al. (1996) and Burrows & Thompson (2002). 

The number emission rate and energy emission rate for electron capture by protons 
(p + e~ n + z/ e )and those for positron capture by neutrons (n + e + — > p + t 7 e )are given by 

q2 poo 

x(^-Q) 2 e(ge _J fcflT + r (Al) 

poo 

x ^-g) 3 e ,E.-J t .T +1 . (A2) 
G^ f°° 

*( E ° + Q) \(E e+ J /kB T + 1 > ( A3 ) 
G 2 f°° 



*( E ° + Q) \(E e+ J /kB T + 1 > ( A4 ) 

respectively, where Gf = 2.302 x 10~ 22 cm MeV -1 is the Fermi constant, qa — 1-23 is the 
axial vector coupling constant, Q ~ 1.29MeV is the difference of rest mass energy between a 
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neutron and a proton, m e is the electron mass, and \x e is the chemical potential of electrons. 
Here, r\ pn and i] np indicate the inclusion of Fermion blocking effects in the nucleon phase 
spaces and are defined as 

Vpn exp(i] n - rj p ) - 1 ' 

^ n p n n (A6) 

nnv exp(r) p - r] n ) - 1 ' 

respectively. Here, n n and n p are the number density of neutrons and protons, respectively, 
and Tj n and r\ p are the degeneracy parameters neutrons and protons, respectively. These 
quantities, i] pn and i] np , are equal to n p and n n , respectively, in the nondegenerate regime 
(where blocking is not important), but when the nucleons are degenerate, those values be- 
come smaller than n n and n p , respectively. 

The number emission rate and energy emission rate for electron capture by nuclei (A + 
e~ — > A' + is e ) depend on the internal proton number Z and the internal neutron number 
N, and are given by 

R A+e -^A' + u e = ^^g 2 JfN p (Z)N h (N)n A 



f°° 1 
x J dE e E e ^ e -my(E e - Q') \ {E ^ )/kBT + l , (A7) 



G F 2 2 



where 



H^a> + v. = ^^9A^N p (Z)N h (N)n A 

poo 1 

x J qi dE e E e ^ e -my(E e - Q'T e{E ^ e)/kBT + v (A8) 

Q' « ^-/ip + A, (A9) 

fo, Z<20 
A^ P (Z) = Iz - 20, 20 < Z < 28 (A10) 
[ 8, Z > 28, 

6, N <3A 

N h (N) = 34<N<40 (All) 

0, N > 40. 

The quantity A 3MeV is the energy of the neutron l/ 5 / 2 state above the ground state. 



-16- 



The number emission rates and the energy emission rates of v e or v e by pair annihilation 
(e~ + e + — > z/j + z/j) are given by 

p _ 1 (^1 + ^2)^ ape 

K e - +e+ ^ e+l7e - - - ^—^e e -e e+ , [All) 



= 2 36 K^ (e - ee++e - ee+) ' (A13) 

„ _ 1 (Ci + C 2 )^ y - q-qc 



= 2 9 K^ (6e - ee++ee - ee+) ' (A15) 

where the weak interaction constants that appear in the above expressions are 

a = 1-76 x Kr 44 cnr 3 , (A16) 
(Ci + Ca)^ = (CV - C A ) 2 + (C v + C A ) 2 , (A17) 
(Ci + C 2 )^ = (C v - C A f + (C v + C A - 2) 2 , (A18) 

with CV = 1/2 + 2sin 2 w , s\n 2 6 w « 0.23, and C A = 1/2. 

The emission rates of creation of z/ e and v e by plasmon decay (7 — > i/j + ^) are given by 

^ ^ (n ^).(S> ^ +i ^ +2 )- (A20) 

where a* ~ 1/137 is the fine structure constant, and 7 ~ 7oa/ (tt 2 + 3r? 2 )/3 with 70 being 
related to the plasma frequency by 70 = hQo/m e c 2 = 5.565 x 1CT 2 . 

The absorption opacities of neutrinos can be derived from the above expressions of 
emissivity. Using Kirchhoff 's law, the spectral-integrated neutrino optical depths are given 
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by 

Rp+e- ~*n+u a + Re-+A->v c +A' + Re-+e+^v e +u e + R*/^v e +i7 e / . nn\ 

T~n,a,u e — , {I\ZO 

TE,a,„ e = , (A24) 

-Rn+e+^p+;7 e + Rn^p+e-+i7 e + -Re-+e+;/ e +z7 e + -^7— ►i/e+i'e / a r)r\ 

1n+e+^p+i7 e 1n~>p+e- +u e + %- +e+ ->v e +v e + Ql-tVe+Ve ( a oc \ 

l~E,a,u e = , [AZb) 

cu v e ,0 

_ Re- +e+ —n> e +i7 e j~ Rj-^Ve+iJ,, { \07\ 

rB,o,^ = — , (A28) 

where n Uu o and u Vi $ are the neutrino number density and energy density, respectively, that 
are spectral-integrated assuming Fermi-Dirac distribution. 

In calculating the scattering opacities of neutrinos, we consider nucleon scattering, elec- 
tron scattering and heavy nucleus scattering, using Fermi-Dirac distribution function in 
spectral-averaging. First, the nucleon scattering opacities are computed as 

-C n 9 Y ( kBT \ 2 F ^n Ul ) rA?q , 

m N \m e c 2 J F 2+j {rj Vi ) 

where C SjP = [4(Cy — l) 2 + 5a 2 ]/24 for a proton and C S) „ = (1 + 5a 2 )/24 for a neutron with 
a ~ 1.25. For j = one obtains the opacities for neutrino number transport, and for j = 1 
the opacities for neutrino energy transport. The variable Yjvjv means the nucleon number 
fraction in which Fermi blocking effects are taken into account, and here we define as 

Ynn = , M(0/ J N 7 fix, (A30) 

1 + (2/3)max(r7 A r, 0) 

where Y N and r/ N are the number fraction and degeneracy parameter of each nucleon, re- 
spectively. 



The electron scattering opacities are given by 



3oo / T 



8 V m P c 2 



2 



1 + 



4 / 



{Cy + C A f+ l -{C V -C A f 



(A31) 



where n e is the number density of electrons, and the weak interaction constants are Cy = 
1/2 + 2$m 2 6 w for electron types, Cy = — 1/2 + 2sin 2 6 1 ^ for \x and r types, Ca = 1/2 for u e , 
and i7 T , and Ca = —1/2 for u e , and u T . 



The effects of neutrino-nucleus scattering will be important in the outer region of the 
disk, where heavy nuclei are not fully dissociated into nucleons. The opacities are given by 

K ^ A ~ UA lQ\m e c-) F 2+j ( Vl/t y 

where ha is the number density of nuclei with the atomic weight A. 

Finally, we mention about the choice of the degeneracy parameters of neutrinos. In 
calculating neutrino emissivities and opacities according to the formulae shown above, one 
should choose the appropriate values of neutrino degeneracy parameters, r\ Vi . We choose 

(A33) 

exp(-r n>I/ J] , (A34) 
exp(-T nil7 J] , (A35) 

with 

C q = ~vT = Ve + V P - Vn- (A36) 

With this choice, one can reproduce the neutrino number density and energy density at 
chemical equilibrium in the neutrino-thick limit. 



(A32) 



Vu x = 0, 

Vu e = CM 1 - 
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Fig. 1. — Profiles of density (upper- left), temperature (upper-middle), electron fraction 
(upper-right), free nucleon fraction (lower-left), entropy per baryon (lower- middle), and elec- 
tron degeneracy parameter (lower-right) for the accretion flow with M = 0.01 (red), 0.03 
(green), 0.05 (blue), and O.O8M sec -1 (purple). 
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Fig. 4. — The ratios of neutrino cooling flux (F_) to viscous heating rate per unit disk 
surface area (F + ) in NDAF with mass accretion rates of 0.08 (red), 0.2 (green), 0.5 (blue), 1.0 
(purple), and 5.0M© sec -1 (cyan). 
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Fig. 5. — Upper panel: The radial profiles of the optical depth with respect to electron 
neutrinos in NDAF with mass accretion rates of l.O(red), 3.0 (green), 5.0 (blue), 8.0 (purple), 
and 10.0M o sec -1 (cyan). Lower panel: The radial profiles of the optical depth (solid) and 
the scattering optical depth by heavy nuclei (dashed) with respect to electron neutrinos. 
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Fig. 6. — Contours of electron fraction (left-top), a-particle fraction (right-top), and heavy 
nucleus fraction (left-bottom) on (r,M) plane. 
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Fig. 7. — Thermal equilibrium curves on the (E, T) plane and (X, M) plane at radial distances 
of r = 100r 9 {red), 30r g (blue), and 10r g (green). 
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Fig. 8. — Upper panel: Total neutrino luminosity (thick-solid), v e luminosity (thin-solid), 
and v e luminosity (dashed) from the hypercritical accretion disk as the function of mass 
accretion rate. Lower panel:The ratio between total luminosity of neutrinos (thick- solid), 
v e luminosity (thin-solid), and v e luminosity (dashed) from the whole surface of NDAF and 
total rest mass energy flowing into the central black hole per unit time (Mc 2 ) as the function 
of mass accretion rate. 



